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It is well known that planets with short orbital periods (< 10days) are common around 
stars like the Sun!-. Stars expand as they evolve, and thus we expect their close planetary 
companions to be engulfed**. However, this phase has never been directly observed. Here, 
we present the discovery of ZTF SLRN-2020, a short-lived optical outburst in the Galactic 
disk accompanied by bright and long-lived infrared emission. The resulting light curve 
and spectra share striking similarities with those of red novae”® - a class of eruptions now 
confirmed? to arise from mergers of binary stars. Its exceptionally low optical luminosity 
(~ 10% ergs!) and radiated energy (~ 6.5 x 10*! erg) point to the engulfment of a planet 
(of 1 — 10 Jupiter masses) by its Sun-like host star. We estimate the Galactic rate of such 
Sub-luminous Red Novae (SLRNe) to be ~ 0.1—few yr‘. Future Galactic plane surveys are 
well-poised to routinely identify them, revealing the demographics of planetary engulfment 


and the ultimate fate of planets in the inner Solar System. 


Using data from the Zwicky Transient Facility (ZTF) time domain survey!?, we searched 
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for slowly evolving outbursts near the Galactic plane (Methods). We identified a transient optical 
source named ZTF 20aazusyv (hereafter ZTF SLRN-2020) that exhibited a fast rise from quiescence 
to peak outburst flux in ~ 10 days, subsequently fading by ~ 10x over six months (Figure 1). 
The long optical outburst duration together with its faint peak flux distinguishes it from common 
Galactic plane transients resulting from white dwarfs with close binary companions (the ‘dwarf’ 
and ‘classical’ novae!'). The transient also exhibits a mid-IR brightening starting ~ 7 months prior 
to the optical outburst, together with bright mid-infrared (IR) emission (> 50x brighter than the 
optical r-band at ~ 4 months after optical peak) that lasted for = 15 months. No X-ray emission 
was detected in follow-up observations during the outburst using the Swift telescope!?, ruling out 


an unstable disk accretion episode around a neutron star or black hole!’ (Methods). 


The bright mid-IR emission during the outburst is suggestive of emission from a warm 
dust shell surrounding the stellar photosphere. We model the optical to mid-IR spectral energy 
distribution (SED; Methods) at 120 days after the optical peak. The analysis reveal a relatively 
hot inner photosphere (~ 9000 K) surrounded by a warm dust shell with temperature ~ 1000 K, 
located behind a dust visual extinction column of Ay ~ 3.6mag (Figure 2). Using the 90% 
confidence interval on the foreground extinction together with three-dimensional Galactic dust 
distribution maps (Methods), we infer the source to be located at a distance of 2 to 7kpc. A 
joint analysis of the overlap between the different dust maps suggests a likely distance of ~ 4 kpc 
(Methods). Performing the same analysis at ~ 320 days after peak, we find the SED to have 
predominantly shifted into the IR bands caused by an increase in the dust optical depth, that can 


be attributed to the formation of ~ 109M, of dust (for a distance of 4 kpc). 
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Using the best estimate for the foreground extinction, we construct a bolometric luminosity 
light curve for the outburst (Figure 2; Methods). The light curve is characterized by an initial 
plateau at a luminosity of ~ 103? x (d/4 kpc)? erg s”! lasting ~ 25 days (Methods) before fading 
by a factor of 5 over the next ~ 100 days. The effective temperature of the photosphere stays 
constant at ~ (6--7) x 103 K during the plateau phase, presumably regulated by the recombination 
temperature of hydrogen!", before fading and cooling to ~ 5 x 103 K. The total radiated energy is 


~ 6.5 x 1011 (d/4 kpc) erg over the first ~ 150 days. 


On UT 2020 November 20, we obtained an optical spectrum of the transient using the Keck-I 
telescope. The spectrum (Figure 3) exhibits a nearly featureless red continuum containing only 
atomic (Na, Ba, H and Mg) and molecular absorption features (VO and TiO). The lack of such 
atomic emission lines during the outburst is inconsistent with common Galactic plane transients 
such as accretion events in young stars!” as well as accretion/thermonuclear outbursts in dwarf and 


classical novae!® !7 


, where hot gas produces emission lines due to atomic recombination and ions 
(Methods). Instead, the molecular features are suggestive of a cool outer photosphere consistent 
with a M4-III type giant star and effective temperature of ~ 3600 K (Methods). Contemporaneous 
near-infrared (NIR) spectra (Figure 3) obtained with the Palomar 200-inch telescope show only 
broad molecular absorption bands from H20 and likely TiO, VO and CO, consistent with a M7-IlI 
type giant possessing an extended, cool envelope (Methods). Late-time NIR spectra (Figure 3) 


obtained ~ 700 days after the optical peak with the Magellan Baade and Keck-II telescopes show 


only a relatively featureless continuum with broad H20 absorption. 


84 The distinctively cool molecular spectroscopic features together with the substantial reddening 
ss of the SED during the outburst are reminiscent of the class of “red novae”” 13?!, Direct photometric 

ss observations of the binary orbital decay prior to the red nova V1309 Sco’ provide strong evidence 

s7 for these events to be associated with catastrophic mergers of binary stars. As the primary star 

ss engulfs its companion, a powerful outflow is launched from the binary system that gradually cools 

sæ and powers a ‘plateau’ in the light curve of the optical transient via recombination of hydrogen”. 

90 Subsequently, the expanding envelope cools and forms dust leading to the emergence of a photosphere 


«1 dominated by molecular absorption and a luminous, long-lived infrared transient”. 


92 The photometric and spectroscopic properties of ZTF SLRN-2020 share striking similarities 
ss with both Galactic and extragalactic red novae (Methods). Yet, its remarkably low luminosity, 
ə even if placed on the far side of the Galactic disk (< few x 10°° erg s~'), makes it exceptional in the 
ə population of red novae that reach the Eddington luminosity of the primary stars (~ 1033 erg s7! 
æ for a 1M, star). We show ZTFSLRN-2020 in the phase space of luminosity and timescales 
97 (Figure 4; Methods) for red novae together with analytical contours for the mass and velocity 
ss of ejecta in models of stellar mergers®. The extremely low luminosity is reasonably explained 
ə by only ~ 10? — 104M, of hydrogen launched in the outflow. The small ejected mass is 


too consistent with the non-detection of radio molecular line emission in follow-up observations with 


101 the Submillimeter Array (SMA) and Very Large Array (VLA; Methods). 


102 Using a high spatial resolution image obtained with the Gemini-South telescope ~ 2 years 


tos after the outburst peak, we identify a faint progenitor source (Methods) in archival NIR images 


04 


from the United Kingdom Infrared Telescope (UKIRT) Galactic plane survey”. Although limited 
by the photometric errors, its brightness and colors are consistent with a ~ 0.8 — 1.5 Mo star on 
the main sequence or early in the sub-giant branch (radius 1 — 4R,). The infrared progenitor is 
thus similar to the Sun, and to the primary star of V1309 Sco, which was a merger event involving 


a low-mass q = 0.1 companion, where q is the companion to primary mass ratio?”???, 


We draw constraints on the mass of the merging object from both the light curve and pre-outburst 


detections. ZTF SLRN-2020’s ejecta mass and radiated energy are both ~ 10? lower than that of 
V1309 Sco (Figure 4), implying a merger with a very low mass companion of ~ 0.1--1M,. In their 
evolution towards coalescence, merging systems lose mass to their circumbinary environments, 


28.303? We compare the 


mediating the angular momentum loss that drives their orbital decay 
pre-outburst limits and detections to these models (Figure 4). The pre-outburst dust model that 
best matches the data has q = 107? or a companion mass of 10M; (Methods). Therefore all our 


estimates squarely point to a close (< 1d orbital period) substellar companion to the primary star, 


plausibly a Jupiter-like planet strikingly similar to known systems. 


The overall duration of the light curve suggests a mass ejection velocity of roughly 30kms~! 


(Figure 4), which is substantially lower than the stellar escape velocity. We can unify these 
estimates with our theoretical understanding of how planetary engulfment might affect a host star. 
In particular, the smaller the engulfed companion, the less dramatic the disturbance to the primary 
star, and the smaller fraction of material that is expected to be ejected at high enough velocities to 


become unbound***°. The short-lived ~ 25 d plateau phase in ZTF SLRN-2020 may be powered 
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by the ejection and unbinding of a small amount (< 104M) of mass at velocities approaching 
the stellar escape velocity (~ 100kms~!; Figure 4). The radius of ZTF SLRN-2020 (Figure 2) 
remains roughly constant at ~ 3 x 10" x (d/4kpc) cm during the plateau and recedes during 
the fading phase, similar to that expected for the gravitational contraction of a merger remnant”. 
These features suggest that the late-time decay over the next ~ 100d is powered by hydrodynamic 


and thermal readjustment of the star following the ingestion of its planetary companion. 


Our interpretation of ZTF SLRN-2020 as the engulfment event of a planetary mass object 
by a Sun-like star provides evidence for a missing link in our understanding of the evolution 
and final fates of planetary systems. It has been long known that the population of gas giants 
in short orbital periods (‘hot Jupiters’!*:3”) have sufficiently low orbital angular momentum such 
that they are unstable to tidal dissipation and are bound to merge with their host stars**“?. This 


43,44 as well as the 


is consistent with the lack of old planetary systems with short orbital periods 
dearth of close planets around sub-giant stars“. Therefore, the observations reported here offer 
the first direct insight into the effect of planetary engulfment on their host stars to interpret common 
indirect techniques used to infer past planetary engulfment via its effects on the long-term stellar 


6.48. chemical enrichment”?, and stellar rotation??93. With empirical and theoretical 


luminosity 
rate predictions ranging from 0.1 to fewyr"! (Methods) for similar “sub-luminous red novae’ 
(SLRNe), upcoming combined optical and infrared surveys of the Galactic plane may reveal 


many similar events, probing the demographics of planetary engulfment and testing theories of 


the co-evolution of stars and their planetary systems. 
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Figure 1: The discovery location and multi-color light curves of ZTF SLRN-2020. (a) A faint 
progenitor identified in archival NIR images from the UKIRT Galactic Plane Survey, (b) the mid-IR 
transient source detected in NEOWISE images from 2020 , (c) a NIR composite follow-up image 
(the white square masks out a nearby region with a detector artifact) of the transient during 2021 
(Panels (a)-(c) have identical spatial scales) and (d) a zoomed-in high spatial resolution image of 
the IR remnant in 2022. (e) Multi-color light curves of the outburst from the ZTF!°, ATLAS’ 
and NEOWISE® surveys (as indicated in the legend). Upper limits are shown as symbols with 
downward arrows. J and O indicate the times of the NIR and optical spectroscopy of the transient, 
while W shows the epoch of the P200 NIR imaging. The vertical dashed lines show the time ranges 
used to perform SED modeling during the outburst in 2020 (SED Outburst) and after the optical 
transient had faded away in 2021 (SED Late). The inset shows a zoom-in of the early time light 
curve (shaded in grey in the main panel), showing faint i-band precursor emission and detection of 
mid-IR emission prior to the onset of the optical outburst. 
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Spectral energy distribution evolution and bolometric light curve of 


ZTFSLRN-2020. (a) Best-fit model (parameters are shown) for the SED of ZTF SLRN-2020 
(shown as black circles) at ~ 125 days after outburst peak. The black solid lines show the total 
flux, the brown dot-dashed lines show the dust emission, while the green dashed and blue dotted 
lines show the scattered and attenuated stellar emission respectively. (b) Same as (a) but for 
~ 320 days after outburst peak. (c) The bolometric luminosity (top), temperature (middle) and 
radius (bottom) evolution of ZTF SLRN-2020 for an estimated distance of 4 kpc and a foreground 


interstellar extinction of Ay 


3.6mag. The red squares in the top panel show the luminosity 


estimated from the two epochs of DUSTY modeling. For comparison, we also show the evolution 
of these parameters for previous red novae, as well as the L œ t~*/° luminosity decay expected for 
a merger remnant. The radius and luminosity of archival events have been scaled as indicated. 
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Figure 3: Optical and infrared spectra of ZTFSLRN-2020. 
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(a) Optical spectrum of the 


outburst (black) obtained ~ 180 days after peak. The spectrum shows clear atomic and molecular 
absorption features, similar to the M4-III type giant HV 2255. For comparison, late-time optical 
spectra of a previous Galactic red nova (V838 Mon; in orange) and an extragalactic Luminous Red 
Nova (AT 2018bwo; in magenta) are shown after applying the inferred foreground extinction". 
The inset shows a zoom-in of the spectra around the region of the Na D, Ha and Ba II atomic lines. 
(b) The NIR spectrum of ZTF SLRN-2020 (black) at ~ +160 days and ~ +690 days after optical 
peak, showing clear broad molecular absorption features of H20, TiO, VO and likely CO, similar 
to the M7-III giant HD 108849 (shown in red). Similar features are also seen in the NIR spectra of 
the previously known extragalactic red nova AT 2018bwo (shown in magenta). 
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Figure 4: Comparison of the ZTF SLRN-2020 outburst to the phase space of luminosity 
(Loo) and timescales (tgo) for previously known red novae. (a) Contours show inferred ejecta 
physical parameters? — ejecta mass in units of solar and Jupiter mass on the left axis, and outflow 
velocity in units of kms”! shown in magenta on the right axis. The solid and hollow star show 
ZTF SLRN-2020 via its light curve plateau duration and time taken to release 90% of the total 
radiated energy respectively (Methods).The inferred ejecta mass of ZTF SLRN-2020 is a hundred 
times lower than any other known red nova, and the characteristic velocities are also nearly an order 
of magnitude lower. (b) Comparison of the pre-outburst dust mass estimated from the precursor 
mid-IR emission of ZTF SLRN-2020 (shown as stars) to models of pre-coalescence mass loss for a 
1 Mo star evolving off the main-sequence with binary mass ratio q = 103 — 10-1. The solid lines 
correspond to an initial radius of 2 Rọ, while shaded region shows variations for radii of 1 — 4 Rọ. 
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Methods 


Identification in ZTF The Zwicky Transient Facility is a wide-field optical time domain survey 
running out of the 48-inch Schmidt telescope (P48) at Palomar observatory!”. With a field of 
view of 47 square degrees, the instrument achieves a median limiting magnitude of r ~ 20.5 mag 


in 30 s exposures of the sky and a survey speed of ~ 3750 square degrees per hour??? 


. provide 
an overview of the various ZTF surveys undertaken in Phase-I of operations, and the survey 
scheduling system designed to carry out operations to maximize volumetric survey speed. The 
transient source ZTF 20aazusyv (hereafter ZTF SLRN-2020) was first detected in the ZTF difference 
imaging pipeline on UT 2020-05-16 (MJD 58985.37) at J2000 coordinates a =19:09:39.78 and 
ô =+05:35:04.25, and at a magnitude of r = 20.03 + 0.19mag. The corresponding Galactic 
coordinates are | = 39.979890, b = —1.48927, placing this source in the direction of Galactic 
disk. The transient was not detected on MJD 58983.38 to a 50 limiting magnitude of r = 20.48. 


We nominally adopt the mid-point of the last non-detection and the first detection of the transient 


in r-band (MJD 58984.38) as the time of start for the eruption. 


We identified ZTF SLRN-2020 as a candidate classical nova as part of an ongoing systematic 
search for Galactic novae in optical/near-infrared (NIR) time domain surveys®. Nova candidates 
were identified in the public ZTF alert stream (g and r filters) by selecting for slow-evolving, large 
amplitude Galactic plane transients through a custom filter implemented on the kowalski time 


domain astronomy server’. The selection criteria were as follows: 


1. The transient should be located within 10 degrees of the Galactic plane. 
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2. The transient should have no prior history of outbursts in the ZTF alert archive. 


3. The transient should have brightened by > 3 mag from the nearest counterpart in the PanSTARRS1 
(PS1) source catalog, if one exists within 2 arcsec; or the source should have no known 


PS1 counterpart within 2 arcsec. 


4. The transient should exhibit a slow evolving light curve with t2 > 30 days, where tz is the 
time taken to fade by 2 magnitudes from the peak of the outburst. This criteria distinguishes 
classical novae from the majority of dwarf nova outbursts that exhibit fast evolving light 


curves with tə < 30 days!!. 


Photometry from ground-based time domain surveys We retrieved gri photometry of the transient 
from data taken with the P48 ZTF camera!”, that were processed with the ZTF data processing 
system’. Light curves were extracted using forced point spread function (PSF) photometry” at 
the location of the transient in the difference images’', where the location was determined from 
the median of the positions in all published alerts of the transient. We report detections in the 
forced photometry for epochs where the signal-to-noise ratio (SNR) is higher than 30, while 50 
upper limits are reported for other epochs. While the g and r-band data were acquired as part of the 
ongoing ZTF public survey, the i-band data was acquired as part of an internal collaboration survey 
where the reference images were contaminated by the transient emission and therefore non-trivial 
to remove. We thus chose to use i-band photometry reported as part of the PSF-fit source catalogs”? 
for each observation, where the transient is detected starting from the first i-band observation of 


the field on UT 2020-04-16 until the end of the observing season. As the transient was brightest in 
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i-band, we use the peak of the i-band light curve as the reference phase for all other observations. 


Given the high cadence of the ZTF data as well as the relatively slow evolution of the 
transient, we performed an inverse variance weighted binning of the flux measurements in bins of 
3 days to improve the SNR and report them in Table 1; however, we also show the 1-day cadence 
early time ZTF 7-band light curve in Figure 1. We perform the same binning for reporting upper 
limits, where we use inverse variance weighted flux uncertainty to report the 50 upper limit for 
each epoch. We also retrieved forced photometry at the source location from data taken with 
the Asteroid Terrestrial-impact Last Alert System (ATLAS) survey?” using the public photometry 
service’. We performed binning using the same method on the ATLAS data and the transient 
emission is clearly detected in both the c and o filters. The ZTF and ATLAS photometry are 


presented in Table 1 and shown in Figure 1. 


Mid-IR NEOWISE Photometry The transient location was observed during the ongoing NEOWISE 


all-sky mid-IR survey in the W1 (3.4 um) and W2 (4.6 um) bands??? ”3, We retrieved time-resolved 


t576. Due to the location of 


coadded images of the field created as part of the unWISE projec 
the source in a crowded field with nearby blended sources (Figure 1), we used a custom code”” 
based on the ZOGY algorithm”! to perform image subtraction on the NEOWISE images using the 
full-depth coadds of the WISE and NEOWISE mission (obtained during 2010-2014) as reference 
images. Photometric measurements were obtained by performing forced PSF photometry at the 


transient position on the subtracted NEOWISE images until the epoch of the last NEOWISE data 


release (data acquired until December 2021). Transient mid-IR emission is clearly detected in the 


'nttps://fallingstar-data.com/forcedphot / 


23 


362 


363 


364 


365 


366 


367 


368 


369 


370 


371 


372 


373 


374 


375 


376 


377 


378 


379 


380 


381 


subtracted images starting ~ 250 days prior to the optical outburst peak. The template-subtracted 


NEOWISE photometry are provided in Table 2 and shown in Figure 1. 


Swift Observations We obtained X-ray follow-up of the transient with the Swift X-ray telescope 
(XRT’®) on the Neil Gehrels Swift Observatory (Swift!?). The Swift observatory observed the 
location of the transient on UT 2020-11-04 (~ 160 days after i-band peak) for a total exposure 
time of 1.35 ks. No source was detected at the location of the transient down to a 30 limiting flux 
of 8.3 x 103 count s71, corresponding to a 0.3-10 keV flux of 2.8 x 10713 ergs cm” s~! for 
a photon index of I’ = 2. In the same observation, no source was detected with the Ultraviolet 
Optical telescope (UVOT®”) in the UVW2 filter, down to a 30 limit of 22.69 AB mag. The X-ray 
non-detection constrains the X-ray luminosity of the source to be < 5.3 x 1032(d/4 kpc)? erg s7! at 
this phase, and therefore is much fainter (even if placed at ~ 20 kpc) than accretion events around 


neutron stars or black holes that reach at least > 103? ergs~! during the outburst!???. 


NIR imaging follow-up We obtained NIR imaging of the transient location on UT 2021-03-31 
and UT 2022-07-23 with the Wide-Field Infrared Camera (WIRC®?) on the Palomar 200-inch 
telescope. We obtained dithered exposures of the field for a total exposure time of 495 s, 330 s and 
330 s respectively in the J, H and K's filters respectively on each epoch. The data were reduced 
using standard techniques including dark subtraction, flat-fielding, sky subtraction, astrometric and 
photometric calibration”. For the image taken in 2021 March when the transient was bright (Figure 
1), we performed aperture photometry on the reduced images at the position of the transient”. The 


resulting magnitudes are J = 16.64 + 0.02, H = 15.42 + 0.02 and Ks = 14.67 + 0.01 on UT 


>We provide all IR magnitudes in the Vega system. 
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2021-03-31. 


Given the crowded nature of the field as well as marginal evidence for blended sources 
near the transient location in the WIRC images, we obtained additional high spatial resolution 
Adaptive Optics (AO) assisted imaging using the Gemini South Adaptive Optics Imager (GSAOT) 
on the Gemini-S telescope”? 35. The source was observed on UT 2022-04-15 as part of a Director's 
Discretionary Program (GS-2022A-DD-102; PI: K. De). We obtained dithered exposures of the 
field using Laser Guide Star (LGS) correction for a total exposure time of 300s on source and 
300s off the source (for background subtraction) in J-band*. The raw images were detrended 
and stacked using the DRAGONS pipeline”, using the source catalog from the WIRC images for 
astrometric and photometric calibration. The transient is clearly detected in the GSAOI images 


together with an unrelated nearby point source ~ 0.6” away (Figure 1). 


We use the GSAOI images to measure the best position of the transient source to be at 
J2000 coordinates a =19:09:39.783 and ô = +05:35:04.269. The coordinates of the unrelated 
nearby source are a = 19:09:39.815 and 0 = +05:35:04.064. The typical positional uncertainty 
is estimated to be < 0.05”by cross-matching against the Gaia DR3 catalog. We perform aperture 
photometry at the location of the transient to obtain a flux of J = 19.17 + 0.05 mag. We use the 
GSAOI positions to derive flux measurements for the transient source detected in the multi-band 
WIRC images from 2022 July. The source is clearly blended with the unrelated nearby source 


in the WIRC images, and we fit the image around the position of the transient with a tractor 


3Observations were also obtained in the H and K filters with GSAOI but were unusable due to poor atmospheric 


conditions. 
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model?” to obtain deblended flux measurements. We fix the position of the transient and nearby 
source as measured from the GSAOI images, and jointly fit the flux of the two sources using the 
measured PSF of the image to minimize the residual y? using the Markov Chain Monte Carlo 
library emcee®*. The resulting flux measurements for the transient are J = 19.17 + 0.15 mag, 


H = 17.75 + 0.09 mag and Ks = 17.13 + 0.10 mag. 


Archival optical/IR observations The location of the transient was covered prior to the outburst 
in the United Kingdom Infrared Telescope (UKIRT) Galactic plane survey?” between 2007 and 
2011. A faint progenitor is detected at the position of the transient as well as reported in the UKIRT 
source catalog. However, the source is clearly blended with the unrelated nearby source detected 
in the GSAOI images, and the photometric measurements are contaminated with the nearby source 
as confirmed by the shifted position of the flux centroid with respect to the transient position in the 
GSAOI images. In order to obtain improved flux measurements for the progenitor counterpart, we 
fit a tractor model”” to the UKIRT image as in the case of the late-time WIRC images. The 
source is detected in the H and K bands, but fainter than the 30 detection threshold in J-band: 
we tabulate the measured fluxes from the tractor modeling in Table 3. The source location was 
also covered in optical images acquired in the PanSTARRSI survey”, but the progenitor was not 
detected. We obtained 50 upper limits on the progenitor optical flux by querying the PS1 source 
catalog in a 2’ region around the position and estimating the median magnitude of sources at the 


5o detection threshold. The flux limits are reported in Table 3. 


Spectroscopic follow-up Following the identification of the transient, we obtained one epoch 


of optical spectroscopy and one epoch of near-infrared spectroscopy of the outburst using the 
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Low Resolution Imaging Spectrometer (LRIS®’) on the Keck-I telescope and TripleSpec”! on the 
Palomar 200-inch telescope respectively. The LRIS data were reduced using standard techniques 
using an automated pipeline’, while the TripleSpec data were reduced using the spextool 


ipeline”? followed by telluric correction and flux calibration using xtellcor™ 
Pip y g 


. The spectra 
obtained during the outburst are summarized in Table 4. We obtained two epochs of near-infrared 
spectroscopic follow-up of the infrared remnant star ~ 2 years after the outburst using the Folded-port 
Infrared Echellette (FIRE?) on the Magellan Baade telescope and one epoch using the Near-Infrared 
Echellette Spectrometer (NIRES””) on the Keck-II telescope. The observing setups are summarized 

in Table 4. The FIRE data were reduced using the pypeit code” for the echelle mode data and the 

firehose pipeline” for the prism mode data. The NIRES data were reduced as the TripleSpec 

data using the spextool pipeline”, followed by telluric correction using the xt e11cor pipeline”. 
As the late-time data were acquired close in time (Table 4), we performed an inverse variance 
weighted stacking of the spectra to improve the signal-to-noise ratio. The individual spectra were 
flux calibrated to match the contemporaneous late-time JH K photometry obtained from WIRC, 


and the stacked spectrum was binned to the resolution of the Keck-NIRES data (R ~ 2700). The 


final reduced and stacked spectra are shown in Figure 3. 


VLA SiO maser search Two hours of Director’s Discretionary Time was awarded on NSF’s Karl 
G. Jansky Very Large Array (VLA) to observe the field of ZTF SLRN-2020 at 43 GHz using the 
most extended (A-array) configuration on 2022 April 14 under project 22A-464 (PI: De). The 
observations, of which about 53 minutes on ZTF SLRN-2020, were set up to be sensitive to any 


SiO (maser) emission between 42.25 and 43.5 GHz. The SiO transitions were targeted to possibly 
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help confirm the peculiar giant merger product and to possibly measure its line-of-sight velocity to 


get a constraint on a (kinematic) distance to support its remarkable low luminosity. 


The observations used 3C286 to calibrate the flux density scale and the instrumental bandpass 
response. Antenna pointing calibration in the direction of the target field was performed using 
J1851+0035, whereas the gain and reference calibrator J1912+0518 observed 50 cycles of 25 
seconds while interspersing the 65 second observations of ZTF SLRN-2020 using a cycle time 
of 2 minutes (with a slew time of about 10 seconds each way). The data were inspected for 
bad visibilities and further processed using the general PIPEAIPS pipeline reduction procedure 
in AIPS. Calibration and imaging did not reveal any potential issues and the resulting synthesized 
beam measured about 65 by 45 mas at a position angle of —60 degrees. The calibrator J1912+0518 


was measured to have a flux density of about 88 + 2 mJy/beam. 


The setup was in particular targeting the (J-0-1) v=1 and v=2 transitions using a channel 


width of 62.5 kHz (~ 0.5 km/s) to an RMS level of ~ 6 mJy/channel, although all seven possibly 


detectable transitions were covered using 1 MHz (~ 7 km/s) channels with an RMS of ~ 1.2 mJy/channel. 


The ~1.2 GHz line-free continuum bandwidth at 43 GHz was imaged down to an RMS of = 
35uJy/beam. No significant SiO line nor continuum emission was detected in the several arcseconds 
surrounding the infrared position of ZTF SLRN-2020, down to the 3.00 flux uncertainty within 
0.2” of the expected position. For the line emission, we searched up to ~ 60 MHz (+420 km/s) of 


the transition rest frequencies which is a conservative upper limit for objects bound to the Galaxy. 
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Figure 5: Integrated CO line intensity map (2-1 transition at 230.538 GHz; left) and continuum 
intensity map (at 225.538 GHz; right) from Submillimeter Array observations centered at the 
position ZTF SLRN-2020 (shown with a white cross). The white ovals in the lower left corner of 
the images show the shape of the synthesized beam (1.2” x 0.9”; PA = —82°) from the observations. 


SMA search for sub-mm molecular emission 


ZTF SLRN-2020 was observed with the Submillimeter Array on two epochs: 15 March 2022 and 
17 June 2022. The CO 2-1 line and continuum emission maps after combining the data from both 
the epochs, are shown in Figure 5. The SMA was in the extended configuration in both the epochs, 
with 6 antennas operating in the array. The phase center coordinates were RA(2000) = 19:09:38.78, 
DEC(2000) = +05:35:04.3. Both the 230 and 240 GHz receivers were tuned to place the frequency 
of the CO 2-1 line, 230.538 GHz, in the spectral window | in the upper sideband (with each 
sideband providing 12 GHz of frequency coverage). With this tuning, the lower sideband covered 
209.5 to 221.5 GHz and upper sideband: 229.5 to 241.5 GHz. The spectral resolution is 140 kHz, 
but the raw data were smoothed to an effective resolution of 1 MHz per channel. The primary 


gain calibrator was 1830+063, observed periodically for 4 min, following the integration for 15 
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minutes on the target source. The quasars 3C279 and BL Lac were observed for 2 hours each, to 
provide bandpass calibration. Absolute fluxes were calibrated using observations of MWC349a, 
Mars and Callisto. In the first epoch, the system temperature varied from 80 to 200K, with the 
atmospheric zenith optical depth (at 225 GHz) ~ 0.1. The weather was better in the second epoch, 


with 7225 = 0.06 — 0.08, and system temperatures varying from 80 — 150K. 


The SMA data were reduced using the Millimeter Interferometer Reduction (MIR) software 
using the online procedure”. The resulting calibrated visibility data were imaged using the Common 
Astronomy Software Application (CASA)’s tclean task. We imaged a range of channels covering 
a velocity interval of about 100 km/s around the CO 2-1 line, producing a data cube of channel 
maps. The synthesized beam size in the combined observations was 1.25 x 0.85”, with PA of 
104.2°. No source is detected at the postion of the transient. The combined (including both 
datasets) continuum sensitivity is 0.17 mJy/beam while for an individual channel for the CO cube 


is 14.3 mJy/beam over 3 km/s. 


Photometric analysis We compare the high cadence r and 7-band light curves of ZTF SLRN-2020 
to a sample of known red novae from the literature in Figure 6. ZTF SLRN-2020 broadly exhibits 
a fast rise to peak starting ~ 10 days prior to outburst peak followed by a slow decay of ~ 2.5 mag 
over ~ 200 days. While the Galactic objects V1309 Sco and OGLE-BLG-360 exhibit a long and 
slow rise prior to the peak of the light curve, the photometric behavior of ZTF SLRN-2020 around 
peak is similar to V838 Mon. The fast rise and slow decay of ZTFSLRN-2020 is also similar to 


many extragalactic red novae (e.g. AT 2018hso and AT 2017jfs) but the duration is longer than 


*https://lweb.cfa.harvard.edu/rtdc/SMAdata/process/mir/swarm_example/ 
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Figure 6: Comparison of the r and ¿ band light curve of ZTF SLRN-2020 (shown as red and 
yellow circles respectively, as indicated) to visual light curves of Galactic and extragalactic red 
novae from the literature. The light curves have been normalized to peak magnitude. For each 
comparison object, we indicate the photometric band of the archival light curve in the legend. The 
comparison objects include V1309 Sco’, OGLE-BLG-3607!, V838 Mon!?!, M31-LRN-20153!-103, 
AT 2019zhd'™, NGC 3437-OT!®, AT 2017jfs! and AT 2018hso!””, 


the fast decay of AT 2019zhd after peak. ZTF SLRN-2020 also exhibits a rapid decay around 
~ 30 days after light curve peak, similar to AT 2019zhd but with a smaller drop in magnitude. 
Additionally, the pre-outburst i-band variability of ZTF SLRN-2020 is similar to that of V838 Mon 
as well as the nearby red novae in M31 (M31-LRN-2015 and AT 2019zhd) where the progenitor 


variability was detected prior to the main outburst. 


In Figure 7, we also compare the color evolution of ZTF SLRN-2020 to objects in the 
literature. ZTF SLRN-2020 exhibits a rapid reddening around the peak of the outburst followed 
by a slow reddening of the g — r color, while the r — 4 color remains relatively constant. Similar 


reddening of the g — r color is also seen in the literature red novae, although the magnitude of 
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Figure 7: Comparison of the g — r (top panel) and r — 2 (bottom panel) observed color evolution 
of ZTF SLRN-2020 (shown as squares and circles respectively) to a sample of red novae from 
the literature that have contemporaneous multi-band coverage during the outburst. The literature 
objects include some of the objects shown in Figure 6 in addition to AT 2020hat, AT 2020kog!”? 
and AT 2018bwo!'”. In each panel, we show a black arrow indicating the estimated shift in the 
color evolution of ZTF SLRN-2020 accounting for the best estimated line-of-sight extinction. 


the progressive reddening is larger for these objects during the ~ 100 days after peak. The slow 
g — r color evolution of ZTF SLRN-2020 is most similar to AT 2018bwo, which also exhibits 
only a small change in the color over the first ~ 40 days, while the slow r — 2 color evolution 
is most similar to AT 2017jfs. Accounting for the estimated foreground Galactic extinction (see 
Spectral Energy Distribution Evolution), ZTF SLRN-2020 remains relatively blue compared to 


the literature sample. As the progressive reddening with time is attributed to the cooling of an 
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expanding photosphere in red novae followed by dust formation!””, the relatively slow evolution in 


ZTF SLRN-2020 may be related to a lower amount of photospheric expansion and dust formation 


relative to other objects, consistent with the modeling of the bolometric light curve and SED (see 


Spectral Energy Distribution Evolution). 
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Figure 8: Comparison of the optical spectrum of ZTF SLRN-2020 to different types of Galactic 


transients, including a classical nova (PGIR 19brv®, in magenta), a dwarf nova (U Gem 
in yellow) and two EXor type young star 


and a FU-Ori type young star outburst (Gaia 17bi!!! 


110 


in red) 


outbursts (ESO-Ha 99!” in green and ASASSN 15qi!! in blue). The continua of some sources 
have been reddened to match that of ZTF SLRN-2020 for easier visualization. 


512 Spectroscopic features The optical outburst spectrum of ZTF SLRN-2020 (Figure 3) exhibits a 


s13 relatively featureless continuum with only atomic/molecular absorption features of Na, Ba II, Ha, 


54 Mg II, TiO and VO. The TiO molecular bands are characteristic in the spectra of cold, late-type 
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giant stars, with later type giants exhibiting deeper absorption features''*'!>, To identify the 
corresponding spectral type, we compared the spectrum of ZTF SLRN-2020 to a library of stellar 
spectra from the VLT X-shooter spectral library!!?, and found a good match of the TiO/VO line 
depths to the M4-III type giant HV 2255. We applied a foreground extinction of Ay = 5.5 mag” 
to the spectrum of the comparison star using a standard Fitzpatrick extinction law?! and show 
the comparison in Figure 3. The inferred spectral type corresponds to an effective photospheric 


temperature of ~ 3600 K!®. 


The optical spectrum shows a strong absorption line in the NaD doublet. We measure the 
total equivalent width (FW) of the line by fitting a polynomial to the absorption feature, and 
measure the uncertainty by creating 1000 realizations of the spectrum by adding Gaussian noise 
scaled to the root-mean-square noise in the adjacent part of the spectrum. We find EW = 3.65 + 
0.75 A. The strength and profile of the Na D line is known to vary with time in many types of 


explosive transients!!?!?0 


, and has recently been shown to be time-varying in extragalactic red 
novae!?!, The line becomes stronger with time in red novae, likely due to the condensation of dust 
in the envelope. As a result, we are unable to use the Na D feature as an indicator of the foreground 


extinction as for other types of Galactic transients?” !??, 


The spectra of ZTF SLRN-2020 show no signatures of emission lines indicative of hot gas in 
the eruption, as commonly seen in other types of Galactic plane transients. In Figure 8, we compare 


the optical spectrum of ZTF SLRN-2020 to other common types of Galactic stellar outbursts. 


>The Ay used here is consistent with the inferred total dust optical depth modeled in the Spectral Energy 


Distribution Evolution. 
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While classical novae exhibit strong and broad emission lines of the Balmer series and nuclear 
processed material!’, accretion outbursts in binary white dwarf systems (i.e. the dwarf novae) 
exhibit strong Balmer absorption features together with Ha in emission, indicative of the hot gas 
in the accretion process. Similarly, the outbursts of young stars exhibit strong emission lines 
at all phases indicative of hot gas in the outflowing material. Unlike the years-to-decades long 
outbursts of the FU Ori stars’, the light curve of ZTF SLRN-2020 is similar to that of the EXor 


class of young stars that exhibit few month long outbursts!?? 


. However, the optical spectra of 
EXor outbursts are dominated by a forest of emission lines of atomic species like H and Call 
at all phases!*°, indicative of magnetospheric disk emission!?”, unlike that seen in the optical 
spectra of ZTFSLRN-2020 (Figure 8). While red novae also exhibit emission lines at early 
phageg!03:109: 121.128. the lines become progressively weaker with time as the photosphere becomes 


t105 


dominated by molecular absorption features from the newly formed molecules and dust”, similar 


to that seen in the spectrum of ZTF SLRN-2020 around 6 months after outburst peak (Figure 3). 


Figure 3 shows NIR spectra of ZTF SLRN-2020 obtained during the outburst (~ 160 days 
after peak) and after the fading of the infrared transient (~ 690 days after peak). We show a 
zoom-in of the spectra of the individual bands in Figure 9. We identify broad molecular absorption 
features of HO affecting the continua of the H and kK -bands. In addition, we identify bandheads 
of TiO and VO in J-band together with weak !?CO absorption bandheads in K-band. The lack 
of strong emission lines of hydrogen (e.g. in the Paschen and Brackett series) as well as the CO 
band-heads which are ubiquitous in EXor type outbursts!?? !? further distinguishes ZTF SLRN-2020 


from the population of young star outbursts. 
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Comparing the spectra to stars in the IRTF Spectral Library!”’, the broad H2O absorption 
bands and weak molecular absorption features are reasonably matched to a M7-III type star, and 
we show the comparisons in Figures 3 and 9. We note that although the NIR spectrum was 
obtained near the optical spectrum, the NIR spectrum suggests a later spectral type than the optical 
spectrum. Such differences have also been identified in the late-time spectra of previous Galactic 


red novae** !30 


, where the extremely low gravity in the extended envelope of the remnant enables 
the formation of H20 molecules and pushes the condensation of TiO and VO to low temperatures. 
The very late-time spectrum of ZTFSLRN-2020 does not show distinctive absorption features 
seen during the outburst, except for the broad H,O absorption bands in H and K-band. Unlike the 
NIR spectrum in outburst, the late-time spectrum exhibits a clear rising continuum towards redder 


wavelengths in K-band, indicative of warm dust emission as inferred from the modeling of the 


spectral energy distribution. 


Spectral Energy Distribution Evolution The outburst of ZTF SLRN-2020 was detected in multiple 
time domain surveys with wavelength coverage extending from the optical g-band to mid-IR 
WISE-W 2 band. While the optical emission likely arises from a hot photosphere, the bright mid-IR 
emission is indicative of a warm dust shell around the eruption. Here, we model the optical to 
mid-IR SED of the transient to estimate the time evolving properties of the optical photosphere and 


dust shell. We use the radiative dust transfer code DUSTY 1? 13 to fit the multi-wavelength data. 


We assume a spherically symmetric distribution of the dust with a x r~? density profile around the 
star, which is assumed to be a point source. We assume the dust grains to be composed of warm 


silicates as indicated by the O-rich composition of the photospheric spectra!**, and with a MRN 
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Figure 9: Identification of spectroscopic features in the NIR spectra of ZTF SLRN-2020 during and 
after the outburst. In each panel, the gray/black lines represent the raw/binned spectrum during the 
outburst (~ 160 d after peak), while the light brown/brown lines show the raw/binned spectrum 
obtained after the fading of the infrared transient (~ 690 d after peak). We also show a comparison 
with the M7-III type star HD 108849 (in red). Prominent atomic and molecular absorption features 
are marked. 
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grain size distribution! (x a73?) with a minimum and maximum grain size of amin = 0.005 um 
and amaz = 0.25 um. We fix the thickness of the dust shell to be Y = 5x the inner radius (rin) 
of the shell; the model SEDs are found to be relatively insensitive to this assumption since the 


3 — 5 um emission arises primarily from the hotter, inner part of the dust shell. 


We fit the observed SED of ZTF SLRN-2020 at two epochs during the outburst that have 
NEOWISE mid-IR coverage (~ 120 and ~ 320 days after outburst peak; shown in Figure 1) using 


a Markov Chain Monte Carlo (MCMC) wrapper around the DUSTY code!** using the Python 


emcee library®®. We model the foreground wavelength-dependent interstellar extinction using 
a Fitzpatrick law?! extending from the optical to the mid-IR. The resulting free parameters of 
the model are the dust optical depth at 0.55 wm (Ty), the foreground visual extinction (Avy), the 
inner stellar temperature (7!.), the dust temperature at the inner edge of the shell (T4) and the total 
flux (F). We assume flat priors on all the fit parameters and ensure convergence of the posterior 
sampling chains. As multi-color optical detections were available only at the ~ 120d epoch, we 
keep Ay as a free model parameter for this epoch, but fix it at the best derived value from the 
120 d epoch when fitting the 320 d epoch. Nevertheless, we find that the best fit for the 320 d epoch 
is relatively insensitive to the assumed Ay since all the photometric data at this epoch are in the 


infrared bands, where the foreground extinction is less important (e.g. Ax œ~ 0.3 — 0.4 mag). 


The best-fit parameters were derived using the median of the posterior sample distribution 
while their confidence intervals are derived from the 1 6th-84th percentile (68% confidence) interval 
of the distributions. The derived parameters and their uncertainties are listed in Table 5, and shown 


in Figures 10 and 11. During the outburst (~ 120d after peak), the SED is well described by an 
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Figure 10: Corner plots showing the model-fit parameters of the MCMC DUSTY modeling of the 
SED of ZTF SLRN-2020 ~ 120 days after outburst peak. 


inner photosphere with a temperature of ~ 8900K surrounded by a warm dust shell with Ty ~ 
1020K. The optical depth of the dust (ty œ~ 1.8) is relatively low at this epoch. The foreground 


extinction is Ay ~% 3.6 mag, although this parameter is degenerate with the optical depth of the 
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Figure 11: Corner plots showing the fit parameters of the MCMC DUSTY modeling of the SED of 
ZTF SLRN-2020 ~ 320 days after outburst peak (Ay is not used as a free parameter in this fit). 


soo dust shell. While multi-epoch mid-IR data is scarce for the literature sample of red novae, similar 
601 parameters involving a relatively hot photosphere surrounded by a low optical depth dust shell were 


s2 also derived at similar phases in DUSTY modeling for the nearby red nova M31-LRN-2015!. 
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Towards the end of the outburst (~ 320d after peak), the DUSTY fitting suggests that both the 
internal photosphere and the surrounding dust shell have cooled substantially, while the optical 


depth of the dust has increased by a factor of ~ 6. 


We estimate the total mass of the dust shell at an epoch! as, 


Ma _7 { 100 cm? g7! Tin 2 
Pid 3 N06 x 10 1 
Mg i ( pe 1000 R3) V (y) 


where Ma is the dust mass and xy œ% 50 — 100 cm? g7! is the dust opacity for a typical dust-to-gas 
ratio”, In order to derive rin, we assume a distance of 4 kpc, as indicated by the distance estimates 
derived in the Spectral Energy Distribution evolution, and list the corresponding radius and mass 
values in Table 5. The observed increase in the dust shell inner radius between the ~ 120 and 
x2 320 days epochs suggests an expansion velocity of ve; ~ 35kms~'. The low inferred velocity 
of the shell suggests that signatures of such an outflow would be difficult to detect in our low 
resolution spectra. Together with the increase in the optical depth, the increased dust mass suggests 
that dust formation in the ejecta shifts the SED of the transient from the optical to the IR bands. 


-25,103 The mass 


Similar abrupt dust formation has also been observed in previous red novae” 
estimate provides a lower limit to the total (dust + gas) ejecta mass of ~ 10~° Mo, noting that 


only a fraction of the ejecta has likely reached the dust condensation radius within a year after the 


outburst!03: 138, 


Distance constraints and bolometric light curve We use the inferred foreground dust extinction 
to ZTF SLRN-2020 to place constraints on the distance to the source. In Figure 12, we show 


published three dimensional Galactic dust extinction maps at the location of ZTF SLRN-2020 
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Figure 12: Constraints on the distance to ZTF SLRN-2020 using Galactic three dimensional dust 
extinction maps. We show the estimated dust extinction as a function of distance for three different 
extinction maps published in the literature!??-!4', We also show the 90% confidence interval for the 
estimated foreground extinction to ZTF SLRN-2020 based on our SED modeling as the magenta 
shaded region. For each dust extinction map, we show the allowed distance interval within the 
estimated extinction range with shaded vertical bars in the same color. 


(created using the mwdust code'??), the estimated Ay range inferred from the SED modeling 
as well as the allowed distance range for each map within the estimated Ay range. As shown, the 
different extinction maps are roughly consistent with each other within the estimated Ay range. 
Given the possible systematic differences between the different maps, we conservatively derive a 
distance range of ~ 2—7 kpc for the transient, placing this source well within the Galactic disk. As 
the distance ranges suggested by the different maps overlap consistently at a range of ~ 3 — 4 kpc, 


we nominally adopt 4 kpc as the best estimate for the source. 


We use the distance and extinction estimate to derive the bolometric luminosity light curve of 


42 


630 


631 


632 


633 


634 


635 


636 


637 


638 


639 


640 


641 


642 


643 


644 


645 


646 


647 


648 


the optical transient around the peak of the outburst. We use the multi-color photometry from ZTF*, 
binned in epochs of 3-days with coverage in all three (grt) filters, and fit a blackbody function to 
derive the effective temperature, luminosity and radius of the star for an estimated distance of 4 kpc. 
The fitting was performed by y? minimization using the Markov Chain Monte Carlo fitting Python 
library emcee. The results are shown in Figure 2. For comparison, we also show the evolution 
of the corresponding parameters for well studied red novae in the literature” as well as the t~*/° 
luminosity decay expected for the gravitational contraction of an inflated envelope surrounding the 
remnant”. We note that a foreground extinction significantly higher than the estimated interval 


(Ay = 5mag) would imply unphysically high blackbody temperatures (>> 10? K), corroborating 


the estimated extinction and distance range. 


The bolometric luminosity of ZTF SLRN-2020 shows an initial plateau at ~ 103? erg s~! for 
~ 15 days followed by a slow decline to 2 x 10*4 erg s~! over the next ~ 100 days. We estimate the 
duration of the plateau by fitting the bolometric light curve with an analytical model developed for 
light curves of Type II supernovae!* and previously used to model light curves of red novae!™. The 
fitting was performed by x? minimization using the emcee code, and suggests a best-fit plateau 
duration of tp = 25.612? days as measured from the assumed time of eruption. The inferred 


bolometric luminosity from the DUSTY modeling at ~ +120 days after outburst peak is consistent 


with the bolometric fitting at the latest epochs, while the DUSTY modeling at ~ +320 days shows 


evidence of a plateau at very late phases. 


“The ATLAS photometric filters cover very wide bandpasses over the g, r, and 7 bandpasses contemporaneously 


with the ZTF photometry, and hence we do use the c and o-band data for constructing bolometric light curves. 
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We estimate the total energy (E), luminosity (Lgo) and duration (tgo) of the outburst? by 
performing trapezoidal integration on the bolometric light curve. We estimate uncertainties on 
these parameters by creating 1000 realizations of the bolometric luminosity light curve using the 
MCMC uncertainty intervals and then repeating the calculations. The corresponding estimates are 
E ~ (6.541.7) x10"! erg, Log © (8.05.1) x 10*4 erg s~! and too ~ 103-420 days for a distance of 
4 kpc and foreground extinction of Ay = 3.6 mag. Although the tgo duration is found to be similar 
to the plateau duration (tp) in previous objects®, we find tgo ~ 4 x tp likely due to a relatively small 
amount of unbound mass that contributes to the recombination luminosity on the plateau phase 
compared to the late-time gravitational contraction following the plateau. Considering the 68% 
confidence interval for the possible distance and foreground extinction, we find the total radiated 
energy and luminosity to be in the range (1.0 — 12.1) x 10 erg s~! and (1.1 — 13.6) x 10% erg st. 
In particular, even if the source was placed on the farthest side of the Galactic disk (d ~ 20 kpc), 


the estimated extinction would suggest a luminosity of < 3.4 x 1039 ergs~t. 


Progenitor photometry The progenitor of ZTF SLRN-2020 is detected only in the H and K 
filters of archival UKIRT images. Here, we attempt to constrain the progenitor star properties 
using the archival NIR photometry. In Figure 13, we show the position of the progenitor in the 
Mx vs. H - K color magnitude diagram as a function of distance along the Galactic disk using 


different three dimensional dust extinction maps. To constrain the progenitor mass, we also show 


stellar evolutionary tracks for stars with initial masses ranging from 0.8 — 3.0 Mo from the MIST 
database'**. Figure 13 shows that the progenitor colors would be too blue to be consistent with 


any stellar tracks for distances larger than ~ 8 — 10 kpc as per the D03! and M06!” models, 
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Figure 13: The progenitor of ZTF SLRN-2020 in the color magnitude diagram. For different 
distances along the Galactic disk, we use available three dimensional extinction maps to de-redden 
the progenitor photometry fluxes. The results are shown as circles for the D03! map, squares 
for the M06'“° map and as triangles for the G19'4! map. We also show the range of absolute 
magnitudes allowed by the 90% Ay confidence region (from the SED modeling) as blue, red and 
gray shaded regions respectively (see legend). We also plot stellar evolutionary tracks from the 
MIST database for stars of initial masses ranging from 0.8 — 3.0 Mọ. The horizontal bar at the 
bottom shows the estimated uncertainty in the H — K color. 


consistent with our constraints on the distance based on the outburst SED modeling. 


We find that the location of the progenitor within the estimated Ay range (shown as shaded 
region) is consistent with a 0.8 — 1.5 Mo star on or evolving off the main sequence. Within the 


90% confidence interval for the extinction, the photometric colors intersect the 1 Mo stellar track at 
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radii of ~ 1 — 4Ro, which we adopt as the likely initial radius of the progenitor star. However, we 
are unable to further constrain the progenitor properties due to the large error on the photometric 
color. Specifically, the allowed progenitor flux and color are also consistent with a lower mass 
z~ 0.8M,progenitor star; in such a case, the engulfment would have occurred while the star was 


on the main sequence (given the current age of the Galactic disk population) likely driven by 


tidal interactions**“”. The late-time WIRC photometry together with the DUSTY modeling shows 


that the star has both nearly returned to its original photospheric temperature (~ 4000 — 5000 K; 


as suggested by the DUSTY model) while also having faded marginally below the progenitor 


brightness likely due to the optically thick dust shell surrounding the remnant. 


Constraints on pre-outburst dust ZTF SLRN-2020 exhibits a mid-IR brightening in NEOWISE 
data starting ~ 7 months prior to the optical outburst. Here, we use the pre-outburst brightening to 
constrain the evolution of the mass loss rate before the red nova outburst. In Figure 14, we show the 
evolution of the pre-outburst SED of the progenitor. The mid-IR source is clearly detected in both 
W1 and W2 bands = 40d prior to the optical outburst. We derive the dust temperature and mass 
by fitting the W1 and W2 photometry with a single-temperature modified blackbody represented 
by a Planck function multiplied by a grain efficiency factor (including effects of absorption and 


scattering). We use published values for silicate dust absorption and scattering!?”. 


We derive a dust temperature of ~ 650 K and a dust mass of M, ~ (0.26 — 3.20) x 109M, 
over the distance range of ~ 2 — 7 kpc at the ~ —44d epoch. The corresponding SED fits are 
shown in Figure 14. As the transient is only detected in the W2-band for the ~ —244d epoch, 


we estimate the dust mass at this epoch assuming the same dust temperature (650 K; consistent 
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Figure 14: Evolution of the SED of ZTF SLRN-2020 progenitor prior to the onset of the optical 
transient. The orange points show the optical/IR photometry from ~ 6 — 12 years before the 
outburst. The black points show the evolution of the progenitor in the mid-IR starting ~ 1.7 years 
before the transient. Solid points indicate detections while hollow points denote 30 upper limits. 
For the -30d and -244d epochs, we show the best fit Silicate dust emission model fit using the 
mid-IR photometry (see text), while we show estimated upper limits to the dust emission using 
the W2 photometry. For the last epoch, we also show the pre-outburst optical brightening of the 
progenitor, coincident with the mid-IR source detected in NEOWISE. 


6s With the W1 non-detection of the source) and use the W2 photometry to obtain a dust mass of 
es Mg ~ (0.7 — 8.8) x 1077 Mg. Instead, if we assume a lower dust temperature of 400 K, the 


s97 corresponding mass estimate is (1.6 - 19) x 1076 Mo that is unlikely as it is higher than the —44 d 
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epoch; we therefore nominally assume a temperature of 650 K. For prior epochs with no W1 or W2 
detection, we only derive the maximum dust mass consistent with the NEOWISE non-detections; 
the SEDs fit to the NEOWISE non-detections suggest that the UKIRT NIR progenitor photometry 
(= 10 years prior to outburst) is not affected by the mid-IR dust emission. The temporal evolution 


of dust mass is shown in Figure 4. 


The transient also exhibits an optical brightening in i-band (compared to the archival PS1 
limits), detected starting ~ 40d before the onset of the optical outburst. Figure 14 shows that 
the progenitor had brightened by > 20x in i-band since the non-detection in PS1, however, it 
cannot be explained as part of the thermal emission from the dust and likely arises from the stellar 
photosphere. The lack of multi-color photometry precludes constraints on the temperature and 
radius of the star during this phase. Similar behavior, arising out of expansion and cooling of the 
photosphere, or tidal heating prior to red nova outbursts has been suggested in V1309 Sco*®*°, 
M101-20150T!8 and M31-LRN-20153!!9. Although we do not have extensive pre-outburst 
coverage in 7-band where the brightening was detected, the plateau in optical luminosity prior 
to the outburst is strikingly similar to V838 Mon, M31-LRN-2015 and AT 2019zhd (Figure 6), 
where the constant luminosity has been attributed to ejected mass that is partially bound and forms 


a constant radius photosphere around the binary?023:193, 


Observed rate estimates We estimate the observed rate of optical transients like ZTF SLRN-2020 
by performing Monte Carlo simulations of the ZTF observing schedule together with a simulated 
population of ZTF SLRN-2020-like events in the Galactic plane. Our method draws from previous 


work on estimating rates of optical transients in ZTF !*4 together with realistic methods for estimating 
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the space distribution of Galactic plane outbursts®. We begin by constructing a template light 
curve for ZTF SLRN-2020 by fitting a Gaussian process model with a constant kernel to the 
peak-magnitude normalized light curve of ZTF SLRN-2020 in the r-band between phases of —20 d 
and 150 d from light curve peak. We do not include upper limits in the fitting process. We choose 
to perform this analysis only in r-band since the public survey data (where this transient was 
identified) includes only g and r band, while the fainter g-band light curve is poorly sampled at 


phases = 15 days away from peak. The resulting template light curve is shown in Figure 15. 
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Figure 15: r-band template light curve of ZTF SLRN-2020, normalized to the peak magnitude of 
the outburst. The dots with error bars show the observed light curve while the red solid line and 


shaded region shows the best-fit template and its 50 confidence interval. 


Using the normalized light curve template, we simulate the observing schedule of the ZTF 


survey between the period of 2019-06-01 and 2021-11-30 (the period over which this search was 
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performed). We assume that the probability of ZTF SLRN-2020-like outbursts follows the stellar 
mass density profile of the Milky Way and therefore simulate a population of transients following 
the density profile in ! (their Table 1 and 2 for a contracted halo model). For each transient, 
we assume the outburst light curve follows the shape of the normalized template scaled to a 
peak absolute magnitude of W, = 2.0 (as inferred for ZTF SLRN-2020) with a nominal scatter 
following a Gaussian distribution with standard deviation o = 0.5 mag. Using the Galactic spatial 
location of each simulated transient, we estimate the apparent magnitude evolution as a function 


of light curve phase using three-dimensional dust distribution maps !?!. 


As in previous works”, for each value of the Galactic event rate, we simulate a population 
of outbursts such that the number of events follows a Poisson distribution with a mean equal to the 
assumed rate. We perform 100 simulations for each value of the assumed rate, repeating the actual 
ZTF observation schedule to estimate the number of recovered events. Using the same selection 
criteria for long-lived Galactic plane outbursts as used to identify ZTF SLRN-2020, we calculate 
the median number of events that pass our selection criteria within this simulated sample, as well 
as its 68% confidence interval. Figure 16 shows the number of events expected to pass our filter as 
a function of the input Galactic rate. In order to derive the best-fit rate and its uncertainty given the 
one confirmed event in our search, we construct a distribution of the fraction of simulations that 
produce the one observed event as a function of the global rate. We fit a skewed normal function 


to this distribution to estimate a Galactic rate (with 68% confidence intervals) 
ro = 4.3153 yr! (2) 


where ro is the Galactic event rate. We show the best-fit distribution in Figure 16. However, we 
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Figure 16: (Left) Number of simulated ZTF SLRN-2020-like events passing our selection criteria 
for an input Galactic rate of events and using the ZTF observing schedule during the search 
period. Black circles show the mean number of novae detected for an input rate while the shaded 
region corresponds to the measured 16-84th percentile range. (Right) Histogram of the fraction of 
simulations producing the number of confirmed events (= 1) as a function of the input Galactic 
rate. The distribution is fit with a skewed Gaussian function (shown as black solid line) and the 


dashed black line shows the best-fit peak of the distribution. 


caution that i) the one observed event represents a lower limit to the actual number of events 
since it was not possible to obtain follow-up spectroscopy for all transients, ii) rate estimates 
from a single observed event are subject to large uncertainties and iii) the observed rates are 
subject to the recovery efficiency of the ZTF subtraction pipeline for faint transients in dense 
Galactic plane fields which have not been quantified for this simulation. Nevertheless, the observed 
rate estimate provides an order-of-magnitude estimate of the Galactic rate for comparison with 
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theoretical estimates. 


Light Curve Modeling and Interpretation ZTF SLRN-2020 has a uniquely low energy among 
the population of red novae transients. Here we discuss the interpertation of this light curve, with 
a focus on estimating the mass of the unseen companion object whose engulfment caused the red 
nova outburst. It is useful to compare ZTF SLRN-2020 to the properties of the particularly-well 
studied source V1309 Sco. Because of the extended period of pre-outburst eclipsing behavior in 
V1309 Sco, the system properties pre-merger are well known. The primary star of V1309 Sco 
was an ~ 1.4Mo, 3.5Ro sub giant”, and the mass ratio was 0.094”, implying a companion mass 
of ~ 0.13M,. The radiated energy was ~ 3 x 10“ erg. The outburst itself ejected an inferred 
~ 3 x1072M, of material with an observed velocity of 150 km s71, about 1/3 the surface escape 


velocity of the donor star 3. 


By comparison, ZTF SLRN-2020 has very similar inferred primary star properties of ~ 1 Mo 
and a likely radius of 1 — 4Rọ. This association implies a similar primary star at a similar 
evolutionary state. The radiated energy of ~ 6.5 x 1011(d/4 kpc)? erg represents a significantly less 
energetic coalescence than V1309 Sco. Similarly, the inferred ejecta mass is much lower than that 
of V1309 Sco at ~ 3x10~° Mo at 100 kms“! during the plateau phase of the recombination-powered 
transient (Figure 4). In merger-powered transients, the loss of orbital energy powers the outburst. 
Thus, a characteristic energy scale is the gravitational binding energy of the pair separated by 
the donor star radius, AE,:,, ~ GMjM:/R,. Given that Mı and R; are similar for V1309 Sco 
and ZTF SLRN-2020, we can use the observed energetics to trace the companion mass M,. The 


radiated energy of ZTF SLRN-2020 is ~ 10? — 103x lower than V1309 Sco. Similarly the ejecta 
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mass (which is directly proportional to the events’ kinetic energies given the similar system escape 
velocities) is ~ 103x lower in ZTF SLRN-2020 than that of V1309 Sco. Scaling down from 
the companion mass of ~ 0.13Mo, these properties point to a planetary mass companion for 


ZTF SLRN-2020 with mass of 0.1 to 1 Jupiter-masses (0.1 — 1M3). 


There are theoretical suggestions that the lower end of this mass range may be unrealistic. 
In particular, models of planet engulfment suggest that low-mass planets engulfed by their host 
stars may fall in so slowly that they do not produce appreciable disturbances or outbursts’. This 
indicates that the outburst properties might not scale linearly with companion mass down to a 
mass as low as 0.1M73. Additionally, the association of this object with other, merger-powered 
red novae transients strongly suggests this scaled-down but otherwise identical physical process 
of engulfment of the companion. Among gaseous planets, lower-mass objects have lower density, 
such that a planet smaller than ~ 1M, would be disrupted by tides outside of a solar-mass main 


4,160 


sequence star rather than plunging in to shock heat and eject stellar atmosphere material. 


Constraints from Pre-Outburst Detections We can also draw information about the engulfment 
that occurred from the IR detections and limits in the years prior to the coalescence. We base these 
constraints on models of circumbinary mass loss in systems trending toward merger. These models 
have successfully explained the dynamics of orbital decay and the observations of increasing 
obscuration in V1309 Sco?330:32192 As the binary system approaches merger, mass loss near 
the outer, L, Lagrange point progressively drains angular momentum of the binary, causing it to 
tighten and further enhancing the mass transfer rate. We estimate the cumulative mass loss to 


the circumbinary environment based on the orbital change from the Roche limit separation, where 
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mass exchange begins, to the point of engulfment, when the separation equals the donor radius. We 
use the RLOF software!®, which integrates the cumulative mass loss given coefficients of specific 


angular momentum loss motivated by hydrodynamic simulations of this process???. 


To do so, we must assume several parameters of the binary system. We have some constraints 
from the progenitor photometry to guide our choices. A nearly sun-like donor star (approximately 
1M, and 1 to 4R,,) appears to be most consistent with the existing constraints. We adopt a 1M. 
and 2Rọ fiducial model, and compare to 1M,, 1Ro and 1M,, 4Ro examples to demonstrate 
the degree to which varying this within the uncertainty affects the result. We further assume that 
the structure of the star is approximated by a F, = 4/3 polytrope, that the gas adiabatic index is 
y = 5/3 and that the star is not initially corotating with the companion orbit. These parameters 
each moderately affect the specific angular momentum loss that accompanies mass loss and orbital 
decay (~ 50% differences in the predicted mass loss with varying parameters”). We then convert 


the estimated total mass lost to a dust mass assuming that the dust-to-gas mass ratio!” is 5 x 107°. 


A caveat is the extrapolation of these models to several orders of magnitude lower companion 
mass than the binary coalescences in which they have been tested. We caution, therefore, that if 
pre-outburst mass loss were caused by a different mechanism - like tidal heating of the stellar 
surface layers during increasingly rapid orbital decay — we might derive different estimates. These 
dust masses expected from modeled systems are shown in Figure 4. We see that the modeled 
dust masses and non-detections in the pre-outburst data are broadly consistent with a model that 
exhibits an exponentially increasing circumbinary distribution of gas and dust as it progresses 


toward merger. In overall normalization, our curves are most consistent with a mass ratio of q ~ 


54 


817 


818 


819 


820 


821 


822 


823 


824 


825 


826 


827 


828 


829 


830 


831 


832 


833 


834 


835 


836 


837 


10:3 to 10~, implying a companion mass of ~ 10-2Mo, or 10M,. This mass is larger than the 
mass inferred from the recombination powered transient’s mass and energy budget, as discussed 
above. However, both mass estimates are strongly suggestive of a planetary mass companion object 


producing the ZTF SLRN-2020 transient. 


Understanding a Planetary Engulfment The planetary mass range inferred from various lines of 
evidence points to the fact that ZTF SLRN-2020 probes new parameter space among the red novae. 
Hydrodynamic simulations of stellar coalescence have suggested that as the mass ratio decreases, 
a smaller fraction of the shock-heated material is truly unbound from the system and ejected 35-36, 
The remaining, bound material settles into a shock-heated, rotating envelope in which spiral shocks 
play an important role in redistributing energy and angular momentum”. An extension of this is 
the limit that the smallest companions should be too insignificant to eject any mass from the host 
star when they are engulfed. Thus, we suggest that the combination of pre-cursor and lightcurve 
data shown in Figure 4 are indicative of a system in which the majority of the shock-heated stellar 
atmosphere remains bound. Indeed, this is strongly suggested by inferred velocities in Figure 4, 


which are much less than the local escape velocity at 3 x 10!! cm of ~ 300 km st. 


If most of the stellar atmosphere mass remains bound following the engulfment, this would 
attribute the plateau luminosity to the recombination of a comparatively small unbound mass 
fraction (~ 3 x 10?M,). It also informs our understanding of the brightness and colors of 
ZTF SLRN-2020 over the hundred-day observing window. ZTF SLRN-2020’s temperature remains 
higher than other red novae and its photosphere radius, following its initial expansion to ~ 3 x 


10 cm & 4.3Ro, gradually decreases rather than growing. These trends both indicate that 
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s38 following an initial period of mass ejection, we may be witnessing the hydrodynamic and thermal 
s3 relaxation of still-bound stellar material that has been disturbed by the engulfment of the substellar 


-4/5 


sax companion. Indeed, the late-time light curve follows a roughly t power-law, which is consistent 


sı with observations of other merger remnants”. 


842 We note that the small ejected mass is consistent with the non-detection of molecular line 
sas emission in the SMA and VLA observations. While a detailed molecular excitation analysis is 
s4 beyond the scope of the work given the non-detection, we can interpret these data with comparisons 


eas to the molecular line observations of V1309 Sco! 


. Observations with the Atacama Large Millimeter 
sas. Array detected the 230 GHz CO transition in V1309 Sco ~ 8 years after merger at a flux level of 
847 3 3.3Jykms7! (peak line flux of ~ 13 mJy for the measured line width of o ~ 100 km s71). 
s As the distance to V1309 Sco’ is similar to that estimated for ZTF SLRN-2020, and assuming 
sag similar molecular excitation conditions as in V1309 Sco, we estimate that the relevant line flux in 
so ZTF SLRN-2020 would be ~ 103x V1309 Sco (peak line flux of < 0.02 mJy) if the molecular 
s51 line fluxes scale approximately linearly with the inferred mass of the unbound material (Figure 
s2 4). Noting that ZTF SLRN-2020 was observed at a much earlier phase and therefore may have 
a3 different excitation conditions, the non-detection of molecular emission is consistent with the small 
s54 inferred ejecta mass in the eruption. While SiO maser emission at 43 GHz was only detected 
sss for V838 Mon!??!98, we use similar scaling arguments and accounting for the larger distance to 


sss V838 Mon, we estimate the maser line flux in ZTF SLRN-2020 to be ~ 1074x fainter (< 0.5 mJy 


ss7 at peak for a resolution of ~ 0.5 kms” !) and hence consistent with the non-detection (< 6 mJy). 


858 The upper end of the range of mass possibilities for the companion of ZTF SLRN-2020 
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nears the mass range of brown dwarfs, which extends from roughly 13M, to 80M ,. Could 
ZTF SLRN-2020's companion object be a brown dwarf rather than a giant planet? We suggest 
that this is substantially less likely because brown dwarf companions at similar separations to hot 
Jupiters are nearly an order of magnitude more rare”. The close-companion portion of this “brown 
dwarf desert” may be caused by the more-rapid tidal decay of these objects orbits into their host 
stars: even if these objects were formed at similar rates to hot Jupiters, many would coalesce with 
their host stars during the pre-main sequence. Since ZTF SLRN-2020 shows no evidence of being 


a young star there seems to be substantially more support for a planetary companion. 


Theoretical Planetary Engulfment Event Rates Planets and substellar objects decay toward 
merger with their host stars through a combination of tidal orbital decay (due to asynchronous 
rotation of the host star with the companion orbit) and stellar evolution (which increases the radius 
of the host star while lowering its spin rate)”. The currently observed population of hot Jupiter 
exoplanets all have inferred tidal decay times of less than 10!° yr. Tidal decay has been directly 
detected for WASP 12b4!:47:!7°171_ Even without tidal decay, each of these planets will surely 


merge with its host star as these stars increase in radius on the sub giant and giant branches. 


There are uncertainties in the propogation of this observed distribution to an engulfment 
event rate. As Metzger et al* have noted, the distribution of known planets at the closest separations 
implies a much higher coalescence rate than the apparent number at larger separations (see their 


Figure 1). This points to a need to resupply planets to these close-in orbits to maintain the 


currently-observed distribution. The occurrence of very close, Jupiter-like planets implies a coalescence 


rate of ~ 0.1 yr~! in the galaxy“, while the number of hot Jupiters at slightly larger separations 
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(with tidal decay timescales of 10° to 101°) appears to be fewer. Whether this relates to observational 
y pp 


selection effects or is an intrinsic property of the population is currently under discussion”! 


With these caveats in mind, we can produce an order of magnitude estimate of the giant-planet 
engulfment rate as follows. The total mass of stars in the Milky Way is ~ 6.4 x 10?M,!?. Ina 
Salpeter IMF, approximately 1/3 of this mass comes from stars greater than a solar mass (which 
is roughly the fraction that evolve off of the main sequence in less than 10! yr). This implies 
that roughly 10!? stars have evolved off of the main sequence in the history of the galaxy. If 1% 
of these stars hosted hot Jupiters? !”, then the time-averaged event rate would be (10°/10!°yr) ~ 
107? yr. In reality, this time-averaged estimate represents a lower limit because it does not 
account for the apparent resupply of planets to close in orbits, perhaps through dynamical scatterings 
(for example, Stephan et al!” estimate an engulfment rate of 0.1 yr~! from A-type stars alone). The 
normalization of these theoretical rates suggests that the detection of a giant planet engulfment by 
ZTF is perhaps fortunate, but not extremely unlikely. On the other hand, observational constraints 
on this rate through transient discoveries may soon provide important insights into the architecture 


and evolution of extrasolar planetary systems. 


A Missing Link in Star-Planet Coevolution As the existence of exoplanetary systems have 
become increasingly clear’, the discovery of a planetary engulfment event has been long-awaited. 
There have been numerous theoretical predictions of the multi-wavelength signatures of these 
events as they occur* 534.178.179. In particular, predictions have focused on their possible appearance 
as optical transients, but many have over-estimated the luminosity and underestimated the event 


duration relative to ZTF SLRN-2020, because they have imagined fast ejecta moving at similar to 
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the stellar escape velocity. A notable exception comes from Yamazaki et al!” 


, whose predictions 
for emission from a shock-heated plasma bubble show characteristic temperatures, photosphere 
radii, and peak luminosities very similar to those observed in ZTF SLRN-2020, strengthening the 
interpretation of a small quantity of shock-heated mass powering the early emission. More detailed 
comparison of this range of models to ZTF SLRN-2020’s emerging properties will be crucial to 
guide future search efforts. In particular, ZTF SLRN-2020’s high IR luminosity (relative to the 


optical bands) suggests that searching for similar events in the IR bands may be much more fruitful 


than in the optical. 


An even larger body of effort has been devoted to understanding the long-lasting effects of 
substellar engulfment???3, with the motivation that if a large fraction of stars engulf one or more 
planets as they evolve, perhaps these events leave long-lasting impacts on the observable stellar 


characteristics. One possible property that has been invoked is secular declines in luminosity, 


which might trace the cooling of a star following the mechanical addition of heat from an engulfment? “. 


Another is enrichment of the stellar atmosphere by lithium and other elements carried by the planet 
bulk composition in higher abundance than the star”. Lithium is thought to be a particularly useful 
tracer because its fragile structure is dissociated in stars, but persists in substellar objects® 50-52, 
Observations of the remnant of ZTF SLRN-2020 at relatively high spectral resolution could search 
for signs of lithium enrichment and serve as a useful benchmark for these theories, as well as 


compare to recent studies of other red novae remnants!*°, 


A final possible signature that has 
been extensively considered is stellar rotation, with enhancements in otherwise-slow rotation rates 


thought to be be possible due to the deposition of angular momentum as the planet is engulfed??-””. 
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Since the total angular momentum budget of a substellar engulfment is well-known, the impact 
on the observable spin rate (as well as the evolution of that spin rate) becomes deeply revealing 
about the transport of angular momentum in the stellar interior. As future work traces each of 
these properties in ZTF SLRN-2020 and similar future transients, we anticipate that these events 
will serve as a crucial missing link in connecting the properties of observed planetary systems to 


the transients they produce and their effects on their host stars. 
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Table 1: Optical/NIR photometry of ZTF SLRN-2020 in the AB magnitude system. 


MJD Phase (d) Filter Mag Instrument 
58955.44 —38.95 r > 19.70 P48+ZTF 
58972.44 —21.95 r > 20.41 P48+ZTF 
58975.46 -18.93 r > 20.14 P48+ZTF 
58978.49 —15.90 r > 20.25 P48+ZTF 


58983.38 -11.01 r 20.56 +0.18 P48+ZTF 


58986.46 —7.93 r 19.37 +0.05 P48+ZTF 


58991.48 —2.91 r 18.07 +0.03 P48+ZTF 


58995.46 1.07 r 17.95 +0.01 P48+ZTF 


58998.46 4.07 r 18.01 +0.02 P48+ZTF 


59005.31 10.92 r 17.90 +0.14 P48+ZTF 


59009.40 15.01 r 18.41 +0.03 P48+ZTF 


59013.35 18.96 r 18.45 +0.02 P48+ZTF 


59017.39 23.00 r 18.64 +0.03 P48+ZTF 


59022.27 27.88 r 18.96 +0.04 P48+ZTF 


59025.37 30.98 r 18.74 +0.04 P48+ZTF 


59033.34 38.95 r 18.89 +0.07 P48+ZTF 


59037.37 42.98 r 19.02 +0.04 P48+ZTF 


59041.29 46.90 r 19.10 +0.04 P48+ZTF 


59045.34 50.95 r 19.10 +0.04 P48+ZTF 
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Table 1: Photometry of ZTF SLRN-2020 in the optical and NIR bands (continued). 


MJD Phase (d) Filter Mag Instrument 
59049.34 54.95 r 19.07 +0.04 P48+ZTF 
59052.36 57.97 r 19.15 +0.03 P48+ZTF 
59057.33 62.94 r 19.15 +0.05 P48+ZTF 
59061.25 66.86 r 19.06 +0.10 P48+ZTF 
59064.32 69.93 r 19.26 +0.06 P48+ZTF 
59070.30 75.91 r 19.41 +0.05 P48+ZTF 
59075.25 80.86 r 19.37 +0.08 P48+ZTF 
59079.21 84.82 r 19.51 +0.05 P48+ZTF 
59082.25 87.86 r 19.55 +0.10 P48+ZTF 
59087.25 92.86 r 19.42 +0.12 P48+ZTF 
59091.19 96.80 r 19.58 +0.10 P48+ZTF 
59094.30 99.91 r 19.31 +0.08 P48+ZTF 
59098.17 103.78 r 19.41 +0.07 P48+ZTF 
59107.26 112.87 r 19.57 +0.15 P48+ZTF 
59111.21 116.82 r 19.68 +0.05 P48+ZTF 
59114.28 119.89 r 19.64 +0.10 P48+ZTF 
59118.23 123.84 r 19.70 +0.13. P48+ZTF 
59122.20 127.81 r 19.55 +0.11 P48+ZTF 
59129.11 134.72 r 19.91 +0.16 P48+ZTF 
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Table 1: Photometry of ZTF SLRN-2020 in the optical and NIR bands (continued). 


MJD Phase (d) Filter Mag Instrument 
59135.14 140.75 r 20.04 +0.11 P484ZTF 
59139.13 144.74 r 20.19+0.14 P48+ZTF 
59146.11 151.72 T 20.26 +0.26 P48+ZTF 
59151.10 156.71 r 20.17 +0.21 P48+ZTF 
59155.08 160.69 r > 18.01 P48+ZTF 
59164.12 169.73 r 19.86 +0.26 P48+ZTF 
59169.14 174.75 r 20.08 +0.19 P48+ZTF 
59182.07 187.68 r 20.14 +0.15 P484ZTF 
59185.08 190.69 r 20.08 +0.12 P484ZTF 
59188.08 193.69 r 20.27 +0.15 P484ZTF 
59194.07 199.68 r > 19.84 P48+ZTF 
59248.56 254.17 r 20.88 +0.28 P48+ZTF 
59251.57 257.18 r 20.62 +0.22 P48+ZTF 
59254.57 260.18 r 20.53 +0.22  P484ZTF 
59264.54 270.15 r 20.76 +0.14 P48+ZTF 
59267.54 273.15 r 21.12 +0.14 P48+ZTF 
59320.48 326.09 r 21.02 +0.29 P48+ZTF 
59373.39 379.00 r 21.25 +0.32  P484ZTF 
58944.48 —49.91 g > 19.96 P48+ZTF 
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Table 1: Photometry of ZTF SLRN-2020 in the optical and NIR bands (continued). 


MJD Phase (d) Filter Mag Instrument 
58955.48 —38.91 g > 20.38 P48+ZTF 
58963.40 —30.99 g > 21.60 P48+ZTF 
58966.46 —27.93 g > 21.38 P484ZTF 
58974.48 —19.91 g > 20.87 P484ZTF 
58978.42 -15.97 g > 20.66 P484ZTF 
58986.41 —7.98 g 21.41 +0.16 P48+ZTF 
58991.40 —2.99 g 19.87 +0.05 P48+ZTF 
58994.43 0.04 g 19.78 +0.04 P48+ZTF 
58998.42 4.03 g 19.83 +0.05 P48+ZTF 
59002.46 8.07 g 19.89 +0.09 P48+ZTF 
59012.33 17.94 g 20.34 +0.07 P48+ZTF 
59016.36 21.97 g 20.72 +0.08  P484ZTF 
59020.39 26.00 g 20.74 +0.09 P48+ZTF 
59024.36 29.97 g 20.65 +0.11 P484ZTF 
59028.36 33.97 g 20.85 +0.18  P484ZTF 
59034.25 39.86 g > 19.92 P48+ZTF 
59037.25 42.86 g > 20.22 P48+ZTF 
59040.27 45.88 g 21.03 +0.13 P484ZTF 
59043.35 48.96 g 20.97 +0.09 P484ZTF 
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Table 1: Photometry of ZTF SLRN-2020 in the optical and NIR bands (continued). 


MJD Phase (d) Filter Mag Instrument 
59048.36 53.97 g 21.11 +0.15 P484ZTF 
59052.29 57.90 g 21.00 +0.10 P48+ZTF 
59055.29 60.90 g 21.29+0.18 P48+ZTF 
59060.38 65.99 g 20.92 +0.19 P484ZTF 
59064.22 69.83 g > 20.34 P48+ZTF 
59067.32 72.93 g 21.52 +0.25 P48+ZTF 
59072.23 77.84 g 21.49+0.21 P484ZTF 
59076.27 81.88 g 21.37 +0.19 P484ZTF 
59080.21 85.82 g 21.11 +0.21 P48+ZTF 
59084.19 89.80 g 21.47 +0.22  P484ZTF 
59091.23 96.84 g > 20.68 P48+ZTF 
59095.21 100.82 g > 20.74 P48+ZTF 
59098.21 103.82 g 21.36 +0.21 P48+ZTF 
59107.17 112.78 g > 20.52 P48+ZTF 
59111.13 116.74 g 21.61 +0.20 P48+ZTF 
59114.21 119.82 g > 21.00 P48+ZTF 
59117.23 122.84 g > 20.32 P48+ZTF 
59121.17 126.78 g > 20.40 P48+ZTF 
59125.17 130.78 g > 20.02 P48+ZTF 
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Table 1: Photometry of ZTF SLRN-2020 in the optical and NIR bands (continued). 


MJD Phase (d) Filter Mag Instrument 
59131.17 136.78 g > 20.26 P48+ZTF 
59135.12 140.73 g > 21.24 P48+ZTF 
59139.11 144.72 g > 21.07 P48+ZTF 
59146.13 151.74 g > 18.03 P48+ZTF 
59149.18 154.79 g > 20.02 P48+ZTF 
59153.13 158.74 g > 20.12 P48+ZTF 
59157.15 162.76 g > 19.88 P48+ZTF 
59167.14 172.75 g > 20.53 P48+ZTF 
59171.13 176.74 g > 20.38 P48+ZTF 


58962.43 --31.96 a 20.65 +0.26 P48+ZTF 


58964.49 —29.90 a 20.41 +0.09 P48+ZTF 


58968.49 --25.90 4 20.48 +0.16 P48+ZTF 


58971.47 -22.92 i 20.53 +0.17 P48+ZTF 


58972.49  --21.90 i 20.57 +0.17 P48+ZTF 


58976.43  --17.96 a 20.13 +0.16 P48+ZTF 


58983.40  --10.99 a 19.56 +0.18 P48+ZTF 


58985.37 —9.02 a 18.38 + 0.03 P48+ZTF 
58991.38 —3.01 a 17.08 +0.02 P48+ZTF 
58994.39 0.00 a 17.06 +0.02 P48+ZTF 
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Table 1: Photometry of ZTF SLRN-2020 in the optical and IR bands (continued). 


MJD Phase (d) Filter Mag Instrument 
58996.41 2.02 a 17.06+0.04 P48+ZTF 
58997.37 2.98 a 17.13 +0.04 P48+ZTF 
59001.35 6.96 a 17.18+0.05 P48+ZTF 
59004.37 9.98 a 17.30 +0.05 P48+ZTF 
59005.39 11.00 a 17.57 +0.11 P48+ZTF 
59008.46 14.07 4 17.54+0.05 P48+ZTF 
59009.34 14.95 a 17.54 +0.03 P48+ZTF 
59012.47 18.08 a 17.53 +0.02 P48+ZTF 
59017.35 22.96 a 17.61 +0.05 P48+ZTF 
59018.33 23.94 a 17.77 +0.03 P48+ZTF 
59022.33 27.94 i 17.87 +0.02 P48+ZTF 
59026.33 31.94 a 17.78 +0.02 P48+ZTF 
59031.31 36.92 i 18.01 +0.05 P48+ZTF 
59032.29 37.90 i 18.02 +0.03 P48+ZTF 
59035.29 40.90 a 18.00 +0.08 P48+ZTF 
59036.35 41.96 a 18.00 +0.02 P48+ZTF 
59040.35 45.96 a 18.08 +0.04 P48+ZTF 
59042.27 47.88 1 18.11 +0.03 P48+ZTF 
59045.32 50.93 a 18.11 +0.04 P48+ZTF 
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1 Data Availability 


2 Code Availability 


The first author will provide python code used to analyze the observations, and any data used to 


generate figures, upon request. 
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MJD Phase (d) Filter Mag Instrument 
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59112.27 117.88 1 18.64+0.07  P484ZTF 
59114.17 119.78 i 18.65 +0.07 P48+ZTF 
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Table 1: Photometry of ZTF SLRN-2020 in the optical and IR bands (continued). 


MJD Phase (d) Filter Mag Instrument 
59116.18 121.79 1 18.65+0.06 P48+ZTF 
59118.17 123.78 1 18.68 +0.03 P484ZTF 
59121.12 126.73 i 18.70 +0.08  P484ZTF 


59121.13 126.74 a 18.70+0.03 P484ZTF 


59125.13 130.74 a 18.79 +0.07 P48+ZTF 


59125.14 130.75 a 18.77 + 0.03 P484ZTF 
59128.19 133.80 a 18.89 +0.08 P484ZTF 
59128.22 133.83 a 18.84+0.05 P484ZTF 
59137.00 142.61 a 19.16 +0.02 P48+ZTF 
58958.58  —35.81 c > 20.53 ATLAS 
58966.52 --27.87 c > 20.29 ATLAS 
58970.57  —23.82 c > 20.01 ATLAS 
58990.48 —3.91 c 19.12 + 0.05 ATLAS 
58994.47 0.08 c 18.81 + 0.05 ATLAS 
58998.47 4.08 c 18.88 + 0.05 ATLAS 
59022.45 28.06 c 20.23 + 0.20 ATLAS 
59026.48 32.09 c 19.75 + 0.10 ATLAS 
59050.41 56.02 c 20.45 + 0.20 ATLAS 
59054.48 60.09 c > 19.88 ATLAS 
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Table 1: Photometry of ZTF SLRN-2020 in the optical and IR bands (continued). 


MJD Phase (d) Filter Mag Instrument 
59071.35 76.96 c 20.55 + 0.28 ATLAS 
59074.41 80.02 c 20.35 + 0.17 ATLAS 
59078.34 83.95 c 20.19 + 0.19 ATLAS 
59082.33 87.94 c 20.52 + 0.25 ATLAS 
59102.30 107.91 c 20.48 + 0.16 ATLAS 
59106.32 111.93 c 20.66 + 0.24 ATLAS 
59110.27 115.88 c 20.38 + 0.24 ATLAS 


59130.31 135.92 c 20.95 + 0.25 ATLAS 


59138.28 143.89 c > 20.32 ATLAS 
59142.22 147.83 c > 20.23 ATLAS 
58944.60  --49.79 o > 19.10 ATLAS 
58949.48  —44.91 0 > 19.42 ATLAS 
58956.59 --37.80 o > 19.86 ATLAS 
58960.62 --33.77 o > 19.93 ATLAS 
58964.54 —29.85 o > 19.43 ATLAS 
58968.53 --25.86 o > 19.96 ATLAS 
58972.51  —21.88 0 > 20.15 ATLAS 
58976.54 --17.85 o > 19.58 ATLAS 
58980.49  —13.90 0 > 19.68 ATLAS 


71 


Table 1: Photometry of ZTF SLRN-2020 in the optical and IR bands (continued). 


MJD Phase (d) Filter Mag Instrument 
58984.58 —9.81 o 19.45 0.06 ATLAS 
58988.47 —5.92 0 18.41 + 0.04 ATLAS 
58992.50 —1.89 0 17.63 + 0.02 ATLAS 
58996.55 2.16 o 17.67 + 0.10 ATLAS 
59000.46 6.07 0 17.79 + 0.02 ATLAS 
59004.41 10.02 o 18.01 + 0.03 ATLAS 
59008.59 14.20 0 18.18 £0.05 ATLAS 
59016.44 22.05 o 18.22 + 0.07 ATLAS 
59020.39 26.00 0 18.45 + 0.05 ATLAS 
59024.34 29.95 o 18.45 0.04 ATLAS 
59028.44 34.05 0 18.54 + 0.04 ATLAS 
59031.47 37.08 o 18.72 + 0.05 ATLAS 
59036.49 42.10 0 18.67 + 0.05 ATLAS 
59040.53 46.14 o 18.91 + 0.06 ATLAS 
59044.43 50.04 0 18.79 + 0.08 ATLAS 
59048.41 54.02 o 18.82 + 0.06 ATLAS 
59052.45 58.06 0 18.77 + 0.06 ATLAS 
59058.33 63.94 o 18.75 + 0.04 ATLAS 
59061.34 66.95 0 18.94 + 0.12 ATLAS 
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Table 1: Photometry of ZTF SLRN-2020 in the optical and IR bands (continued). 


MJD Phase (d) Filter Mag Instrument 
59064.42 70.03 o 18.86 + 0.13 ATLAS 
59068.42 74.03 o 18.96 + 0.06 ATLAS 
59072.44 78.05 O 18.96 + 0.08 ATLAS 
59076.40 82.01 o 18.99 + 0.07 ATLAS 
59080.36 85.97 O 19.25 + 0.07 ATLAS 
59084.33 89.94 o 19.19 + 0.06 ATLAS 
59088.34 93.95 o 19.30 + 0.12 ATLAS 
59092.33 97.94 O 19.27 +0.13 ATLAS 


59096.36 101.97 o 19.02 + 0.06 ATLAS 


59100.30 105.91 0 19.22 + 0.08 ATLAS 


59104.33 109.94 o 19.23 + 0.09 ATLAS 


59108.29 113.90 o 19.36 + 0.10 ATLAS 


59112.28 117.89 0 19.36 + 0.10 ATLAS 
59115.34 120.95 o 19.42 + 0.11 ATLAS 
59120.32 125.93 0 19.19 + 0.14 ATLAS 


59126.25 131.86 o 19.74 + 0.14 ATLAS 


59129.27 134.88 O 19.39 + 0.12 ATLAS 


59136.28 141.89 o 19.71 + 0.18 ATLAS 


59140.23 145.84 0 19.58 + 0.14 ATLAS 
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Table 1: Photometry of ZTF SLRN-2020 in the optical and IR bands (continued). 


MJD Phase (d) Filter Mag Instrument 
59144.21 149.82 0 19.57 20.17 ATLAS 
59148.22 153.83 o > 18.39 ATLAS 


59154.24 159.85 O 19.49 + 0.10 ATLAS 


59158.19 163.80 o 19.53 + 0.17 ATLAS 


59164.23 169.84 O 19.44 + 0.12 ATLAS 


59184.19 189.80 o 19.49 + 0.16 ATLAS 
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Table 2: Template-subtracted NEOWISE mid-IR photometry of ZTF SLRN-2020. 


MJD Phase Filter Flux 


days mJy 


55296.47 —3697.92 W1 0.16+0.11 


55479.51 -3514.88 W1 0.00 0.08 


56940.99 —2053.40 W1  —0.12 +0.13 


57122.03 —1872.36 W1 0.11 + 0.21 


57300.01 —1694.38 W1  —0.02 + 0.14 


57487.55 —1506.84 W1 0.02+0.21 


57658.98 —1335.41 W1  —0.08 + 0.14 


57851.76 —1142.63 W1  —0.09 + 0.19 


58021.99 —972.40 W1 --0.07 +0.12 


58217.44 —776.95 W1 0.02 + 0.22 


58381.60 —612.79 W1  —0.13 +0.13 


58583.15 —411.24 WI 0.02 t 0.22 


58747.25  -24714 W1 0.07 + 0.14 


58949.95 —44.44 W1 0.26 0.22 


59113.93 119.54 W1 2.21 +0.13 


59314.32 319.93 W1 1.03 + 0.23 


59478.41 484.02 W1 0.42 + 0.12 


55296.47 —3697.92 W2 0.03 + 0.09 
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Table 2: Template-subtracted NEOWISE mid-IR photometry of ZTF SLRN-2020 (continued). 


MJD Phase Filter Flux 


days mJy 


55479.51 —3514.88 W2 0.03+0.08 


56940.99 —2053.40 W2  -0.12 + 0.09 


57122.03 -1872.36 W2 0.00+0.10 


57300.01 -1694.38 W2 —0.16+0.10 


57487.55 —1506.84 W2 —0.12+0.11 


57658.98 -1335.41 W2 —0.12 + 0.09 


57851.69 -1142.70 W2  —0.09 +0.11 


58021.99 —972.40 W2  -0.01 +0.10 


58217.44 —776.95 W2 0.00 0.10 


58381.60 —612.79 W2 —0.15 +0.11 


58583.15 —411.24 W2 0.010.13 


58747.25 —247.14 W2 0.130.10 


58949.88 —44.51 W2 0.47 0.12 


59113.93 119.54 W2  2.60240.11 


59314.32 319.93 W2 1.12 + 0.13 


59478.41 484.02 W2 0.40 0.12 
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Table 3: Archival photometry of the progenitor of ZTFSLRN-2020 along with the time of 


observation. Upper limits are reported at 50 confidence. 


Survey MJD Filter Mag 

UKIRT 54340.36 J > 17.50 
54340.37 H 16.66 0.25 
54997.52 K 16.32 +0.18 

PS1 55744.48 g > 22.03 
55806.58 r > 22:02 
55813.56 a > 21.82 
55925.06 ž > 20.97 
56179.10 y > 19.93 


Table 4: Spectroscopic follow-up of ZTF SLRN-2020. The spectra denoted by | were stacked 


together to obtain the final binned late-time NIR spectrum. 


Date Phase Instrument Exposure time Wavelength Range Resolution 
UT d s A 

2020-10-29 +156 P200 + TSpec 1200 10000 — 24500 =~ 2700 
2020-11-20 +178 Keck-I + LRIS 1200 3200 — 10000 =~ 1000 
2022-03-17' +660 Keck-II + NIRES 2400 10000 — 24500 =~ 2700 
2022-04-03' +677 Magellan/Baade + FIRE/Echelle 3600 9000 — 25000 =~ 6000 
2022-05-25' +729 Magellan/Baade + FIRE/Prism 3000 9000 — 25000 z~ 300 - 500 
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Table 5: Derived dust parameters from the multi-epoch DUSTY modeling of ZTF SLRN-2020. The 


radius inner shell and ejecta mass are derived for an estimated distance of 4 kpc. 


Phase F T Ay Ta Ta Tin Mej 

d 10 “ergcm ?s! (A = 0.55 um) mag K K AU Mo 
+120 10,025 5 1:8103 3.6496 8970+598 1014+33 BB x 1078 
+320 11.54179 13.0434 3.6 (fixed) 4300+1$50 41541 53. N02 x 1077 
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